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I. INTRODUCTION 



The calculation of relativistic corrections to energy levels of atomic systems is usually 
accomplished by using the many-electron Dirac-Coulomb (DC) Hamiltonian with possible 
inclusion of the Breit interaction between electrons. However, such a Hamiltonian can 
not be rigorously derived from Quantum Electrodynamic (QED) theory, and thus gives an 
incomplete treatment of relativistic and QED effects. The electron self-energy and vacuum 
polarization can be included in the DC Hamiltonian , though only in an approximate 
way. A different approach, which is justified by quantum field theory, is to start from a 
well adapted one-electron local potential and build many body perturbation theory. This 
approach allows for the consistent inclusion of QED effects as well as a correct treatment of 
;he so called "negative energy states". It is being pursued by Sapirstein and collaborators 
^, but so far no high accuracy results have been achieved for neutral few electron atoms. 
An alternative approach, which is suited for light atoms, relies on expansion of energy levels 
in powers of the fine structure constant 



E{a) = ^(2) + + i^(5) + ^(6) + o{a') 



(1) 



where i?'-"^ is the contribution of order 



m a 



so E''-^^ is the nonrelativistic energy as given 



by the Schrodinger Hamiltonian H'^'^^ = Hq 



Ho = T.(—-—)+T.T.- 



E^^^ is the leading relativistic correction given by the Breit-Pauli Hamiltonian H^^^ 



(2) 



(3) 



where 
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Am? 



X Pa 



a>b b 
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Tab rt ' 



(Ta ■ (Tb () (rafc) + — [6^-3 



4 r^^ 



r 



ab 
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4 r^^ 



2 i^a ■ fab Xpb-ab- fab X Pa) + {^b ' fab X Pb - ^a ■ fab X Pa) 



(4) 
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is the leading QED correction, which includes Bethe logarithms. It has been first 



obtained for hydrogen, for a review see Refs. ly| • A few years later E^^^ was obtained for 
the helium atom 'T'l, see can be easily 

extended to arbitrary light atoms, and recently calculations of E^^^ have been performed for 
lithium ^, Q| and beryllium atoms E^^^ is a higher order relativistic correction and is 
the subject of the present work. It can be expressed as a sum of three terms, 

(^-^O) a>b b 

where H^^^ is an effective Hamiltonian of order ma^. It is well known that the second order 
correction from the Breit-Pauli Hamiltonian is divergent since it contains, for example, the 
Dirac 6 functions. It is less well known that H^^^ also leads to divergent matrix elements, 
and yet less well known that in the sum of both terms these divergences almost cancel out. 
The additional term containing A is the contribution coming from the forward scattering 
three-photon exchange amplitude which cancels the last divergence in electron-electron in- 
teractions, which leads to a finite result. The cancellation of divergences requires at first the 
inclusion of a regulator, a cut-off in the maximum photon momenta, which is allowed to go 
to infinity when all terms are combined together. 

The first derivation of H^^^ was performed for helium fine-structure by Douglas and KroU 



m 



12| . In this case all matrix elements were finite because they considered only the split- 
ting of nPj levels. The numerical evaluation of this splitting has been performed to a high 
degree of precision by Yan and Drake in 



1^ . Since calculations of higher order relativistic 
corrections when singular matrix elements are present are rather complicated, they were 
first studied in detail for positronium, the electron-positron system. The ma^ contribu- 
tion to positronium hyperfine splitting was first obtained (without annihilation terms) by 
Caswell and Lepage in Q], where they introduced a new approach to bound state QED 
namely, nonrelativistic quantum electrodynamics (NRQED). Although their original calcu- 
lations happened to contain some mistakes the idea of NRQED was very fruitful, because it 
simplified enormously the treatment of bound states. Its use has led to significant progress in 
bound state QED, with the calculation of the coniplete thr ee p hoton-exchange contribution 
of order ma^ to positronium energy levels in and |l7j |. It was shown there, that by 

introducing a regulator, either a photon momentum cut-off or dimensional regularization, 
one can derive and calculate all matrix elements in a consistent way. The agreement between 



these calculations and the other purely numerical calculation based on Bethe-Salpeter equa- 
tion 



18l [19^ justify the correctness of the effective Hamiltonian or NRQED approaches. It 
was quickly found, after the positronium exercise, that a similar effective Hamiltonian H^^^ 
can be derived for the Helium atom. Although the derivation of H^^^ for S- and P-states of 
heU.,„ ,s .athe. st..a,gKtfo™..d S. the el„„,.at,o„ of e,eot..o„-e,eot..o„ ..gulant.es a.d 
the calculation of matrix elements is quite involved. For this reason the first results have 



been obtained for triplet states 2^5*1 in 21 1 and 2^P [22], where electron-electron singulari- 
ties are not present, because the wave function vanishes at ri = r2. Within the dimensional 
regularization scheme Korobov and Yelkhovsky [2^ were able to derive a complete set of 
finite operators and calculate their matrix elements for the 1^5*0 ground state of helium. 
None of these results have been confirmed yet. In this work we present a simple derivation 
of effective operators contributing to H^^^ for an arbitrary state of arbitrary light atoms. 
The results obtained agree for the special cases of the ^S, ^Si and ^Pj levels of helium with 



agree lo] 



the former result in [12, |20[. Since we do not explicitly eliminate here electron-electron 
singularities we were not able to verify the result |23] for the ground state of helium. 

Our derivation consists of three steps. The first step is the Foldy-Wouthuysen (FW) 
transformation of a single electron Dirac equation in an electromagnetic field 24], performed 
to the appropriate level of accuracy. The second step is formal. It is the quantization 
of the electrornagnetic field interacting with the atom, using the Feynman integration by 



paths method 



24l |. The third step is the derivation of an effective interaction through the 



perturbative expansion of the equal-time Green function of the total atomic system. 

II. FOLDY-WOUTHUYSEN TRANSFORMATION 

The Foldy-Wouthuysen (FW) transformation is the nonrelativistic expansion of the 
Dirac Hamiltonian in an external electromagnetic field, 

H = d-n + (3m + eA'^, (6) 

where n = p — e A. The FW transformation S [2^ leads to a new Hamiltonian 

HFw = e''iH~tdt)e-"', (7) 

which decouples the upper and lower components of the Dirac wave function up to a specified 
order in the 1/m expansion. Here we calculate FW Hamiltonian up to terms which contribute 



to ma^ to the energy. While it is not clear here which term contributes at which order, 
we postpone this to the next section where this issue become more obvious. Contrary to 
standard textbooks, we use a more convenient single Foldy-Wouthuysen operator which 
can be written as 



S = 



(/3a-7r- -^/3(a-7f)=^+-^[a-7f, e^°-i9t]+y|. (8) 



2m 

where y is an as yet unspecified odd operator {13, Y} — 0, such that \Y, eA^ — idt\ ~ 
[y, (a • tt)^] 0. It will be fixed at the end to cancel all higher order odd terms. The F.W. 
Hamiltonian is expanded in a power series in S 

HFw^Y.n^'^ + ... (9) 

3=0 

where 

n^^^ = J [is,n^^~^^] for J > 1, (10) 

and higher order terms in this expansion, denoted by dots, are neglected. The calculations 

of subsequent commutators is rather tedious. For the reader's convenience we present a 
separate result for each H^^\ 

7^(1) = l^^.^f -JL^a-7r)''-^[d-7r,d-E] + ^[Y,a-7r] 
m 3m^ Aw? 2 m 

-a-n + -^{d-nf-pY--^[{d-7Tf,eA''-idt]--^d-k (11) 
6m^ om"^ Am^ 



^(2) ^ _^(5.^)2^^^.^^4_^_^(^.^)6^^[5.^^.^] 

2m ^ ^ 3m3^ ^ 18 m^ ^ ' 8m^ ^ ' ^ 
1 r,. ^ ie 



[Y,d-n]--—[(d-ny,d-E] 



2 m 24 m- 

1 /3e ^ ^ 

+ 7rA — a[(^ ■ [{(^ ■ 7?)^, -'idt]] - — — ^ (a • 7ra • + a • -Ea • tt) 
2Am^ iDm"* 

(a ■ tt)^ H (a ■ tt)^ [a ■ n,[d ■ 7r,d ■ E]] 



2m^ 3m^ 16 m^ 

■^((a-7f)2a-^ + a-i(a-^)2) (12) 
8 m'' 

Gm-^ 6m^ 96 m^ 



^^'^ = -^(a-ff)^ + /^(a-7f)^ + -^[a-7f,[a-7f,[a-7f,a-E]]] 



' [d-Tf,id-7ffd-E + d-E{d- 7f)'] + [(a • 7f)^ « • i] 



48m4 L ' V / V / J 24m4 



6 m'^ 6 48 m'^ 



+ ^-^((a-nfd-E + a-E(a-nf) (13) 



[a ■ vf, (a ■ vr)^ a ■ E + a ■ E (a ■ vf)^] [(a ■ tt)^, 5 ■ i?] 



192 m4 ^ " ' ' ' ' 96 m 

+ (a -71)5 (14) 

The sum of T^*^*-*, Eq. Q gives a Hamiltonian, which still depends on Y. Following FW 
principle, this operator is now chosen to cancel all the higher order odd terms from this sum, 
namely: 



Y = —^—r id ■ vr)^ — d ■ E -\ [a ■ vr, [a ■ vr, a ■ E]] 

5m4^ ^ 4m2 24m3 ^ ' 

{{d-nfd-E + d-E{d-nf). (17) 
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Y fulfills the initial ansatz, that commutators [F, eA^ — idt] and [Y, {d ■ rr)^] are of higher 
order and thus can be neglected. The resulting FW Hamiltonian is 

(ff-TrV ((T-vr)^ (a ■ vf)*^ ie ^ ^ 
Hfw = + ^—^ - + \— ^ _ or ■ 7r,a ■ E 

~T^^ (a-7fa-E + CT-Ea-7rj- [ff ■ vr, [ct • vf , [a ■ vf, a ■ E]]] 

16 m"* ^ ^ 128 

i e 
16 m 

where we used the commutator identity 



^^A{^-^n^-^,^-E] + [a.7f,a-E]{a-nf], (18) 



[(a■7^)^a■E] = [a ■ vf, [a ■ tt, [a ■ vf, <? ■ E]]] 

+ 2 ((^ ■ 7f)2 [ff ■ 7f, a ■ ^] + [a ■ 7f, ■ ^] {a ■ Tf)^) (19) 

to simplify Hp\y. Moreover, there is some arbitrariness in the operator S, what means that 
Hp\y is not unique. The standard approach which relies on subsequent use of FW- 
transformations differs from this one, by the transformation S with some additional even 
operator. However, all Hpw have to be equivalent at the level of matrix elements between 
the states which satisfy the Schrodinger equation. 



Let us now study the simple case of an external static potential V = eA". The FW 
Hamiltonian with the help of simple commutations takes the form 



This Hamiltonian is equivalent to the one derived previously in |15| , after use of the identity 



which holds for expectation values on stationary Schrodinger states 0. For the exact 
Coulomb potential V = —Za/r, matrix elements of Hdc become singular. Nevertheless, as 
was shown in one can obtain Dirac energy levels up to order m {Z a)^ by regularizing 
the Coulomb potential in an arbitrary way, and all singularities cancel out between the first 
and second order matrix elements. 

Our aim here is to obtain the Hamiltonian for further calculations of m contribution 
to energy levels of an arbitrary light atom. For this one can neglect the vector potential A 
in all the terms having and in the denominator. Moreover, less obviously, one can 
neglect the term with a ■ Aa ■ E and the term. It is because they are of second order 
in electromagnetic fields which additionally contain derivatives, and thus contribute only at 
higher orders. After these simplifications, Hpw takes the form 




(20) 



mpMp',vm = 4mvv)'\<p) 



(21) 




1 



(vr^ — ea ■ B TT^ 



n^ea- B) 




a ■ pa ■ E + a ■ E a ■ p 




(22) 



From this Hamiltonian one builds the many body Lagrangian density 



C = (f)*{i dt - Hfw)(P + ^EM, 



(23) 



where Cem is a Lagrangian of the electromagnetic field, and with the help of perturbation 
theory calculates Green functions. 
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III. THE HIGHER ORDER BREIT-PAULI HAMILTONIAN 



We consider the equal time retarded Green function G = G{{r'^},t'; {ra},t), where by 
{fa} we denote the set of coordinates for all particles of the system. This Green function is 
similar to that used by Shabaev in 2^. In the stationary case considered here, G = G(t' — t). 
The Fourier transform of G in the time variable t' — t can be written as 

which is the definition of the effective Hamiltonian H^s^E). In the nonrelativistic case 
Hes = Hq. All the relativistic and QED corrections resulting from the Lagrangian can be 
represented as 

Hi — ilQ ^ — -^o I2j — tlQ III — tlQ III — -Ho ^ — -"0 



1 1 



(25) 



E-H^-^{E) E-H,siE) 

where ^{E) is the n-particle irreducible contribution. The energy level can be interpreted 
as a pole of G{E) as a function of E. For this it is convenient to consider the matrix element 
of G between the the nonrelativistic wave function corresponding to this energy level. There 
is always such a correspondence, since relativistic and QED effects are small perturbations 
of the system. We follow here a relativistic approach for the electron self-energy presented 
in ^. This matrix element is 

where 

aiE) = miEM) + EmiEMn) (0n|S(E)|0) + . . . (27) 

,„ II' — Ilin 

Having (j{E), the correction to the energy level can be expressed as 
6E = E-Eo = cr{Eo) + a\E^) a{E^) + ... 

= m{Eom + m{Eo) ^ -^{E^m + (0|s'(Eo)|0) m{E^m + . . (28) 

Since the last term in Eq. ()28|1 can be neglected up to order ma^, one can consider only 
S(£'o). In most cases, the explicit dependence of S on state, through i?o, can be eliminated 
by appropriate transformations, with the help of various commutations. The only exception 

8 



is the so called Bethe logarithm, which contributes only to the order ma^. If we consider 
this term separately, the operator S gives an effective Hamiltonian 



eff 



Ho + E = Ho + H'-^^ + H^^^ + H^^^ + ... 



(29) 



from which one calculates corrections to energy levels as in Eq. (jSj). The calculation of S 
follows from Feynman rules for Lagrangian in Eq. ()23p. We will use the photon propagator 
in the Coulomb gauge: 



1 



/i = 1/ = , 



(30) 



and consider separately corrections due to exchange of the Coulomb Gqo and the transverse 
Gij photon. The typical one photon exchange contribution between electrons a and b is: 



+ 



(2 7r)4i 



G,Ak)\{<P 



faik)e 
1 



i k-fa 



1 



Eo-Ho-k^ + ie 
■ fai-k) e-^'-'' 



3l{-k)e 



-i k-ff, 



(31) 



Eq- Ho - k^ + ie 

where is an eigenstate of Hq and is an electromagnetic current operator for particle 
a. One obtains the exact form of j^(/c) from the Lagrangian in Eq. ()23p . and is defined as 
the coefficient which multiplies the polarization vector e'^ in the annihilation part of the 
electromagnetic potential 



A'^(f,t)~e^e 



^J. i k-r—i k'^ t 



dx + h.c. 



The first terms of the nonrelativistic expansion of j° component are 

i ^ 1 ^ 

f{k) = l + - — a-kxp- - — -k'^ + ... 
4m 



and of the / component are 



j[k) = — + —ax k. 
m 2m 



(32) 



(33) 



(34) 



Most of the calculation is performed in the nonretardation approximation, namely one sets 
= in the photon propagator G^v{k) and j{k). The retardation corrections are considered 
separately. Within this approximation and using the symmetrization k^ ^ —k^, the k^ 
integral is 

1 r dk^ [ 1 1 11 



27Ti 



-AE-k'^ + ie -AE + k^ + ie 
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where we have assumed that Ai? is positive, which is correct when is the ground state. 
For excited states, the integration contour is deformed in such a way that all the poles from 
the electron propagator lie on one side, so it is not strictly speaking the Feynman contour. 
However the result of the integration for excited states is the same as in the above, which 
leads to 

imEoM) = -e'J 0ysG,Ak) {<p\fa{k) e''<''^-'^bl{-k) \<l>) . (36) 

The k integral is the Fourier transform of the photon propagator in the nonretardation 
approximation 

O.Url^f^G,.W^±\-X^ . (37) 



One easily recognizes that in the nonrelativistic limit Goo is the Coulomb interaction. How- 
ever this term is already included in Hq, which means that this nonrelativistic Coulomb 
interaction has to be excluded from the perturbative expansion. Next order terms resulting 
from f and Jlead to the Breit Pauli Hamiltonian, Eq. (0)). Below we derive the higher order 
term in the nonrelativistic expansion, namely the m Hamiltonian, which we call here the 
higher order effective Hamiltonian It is expressed as a sum of various contributions 

= (38) 

i=0,9 

which are calculated in the following. 
5Hq is the kinetic energy correction 

= E tI^ • (39) 



SHi is a correction due to the last three terms in Hp]y in Eq. ()22|) . These terms 
involve only A^, so the nonretardation approximation is strictly valid here. This correction 
6Hi includes the Coulomb interaction between the electron and the nucleus, and between 
electrons. So, if we denote by V the nonrelativistic interaction potential 

v'-E-f^ + EEf («) 

a ' o. a>b b ' 

and by £a the static electric field at the position of particle a 

e£.^-V.V = -Za')^+Y.a'f (41) 

a b^a ab 

10 



then 5Hi can be written as 



3 

^ 61 mr ^ 

is a correction to the Coulomb interaction between electrons which comes from the 
^th ^gj^jji \^ Hpwi namely 

^-^{eV ■ E + ea-iExp-px E)) (43) 

If interaction of both electrons is modified by this term, it can be obtained in the nonretar- 
dation approximation Eq. (jHUj) . so one obtains 

47r 1 



a>fe b 

+4 ((?„ X p;) 



Xij 1 / », j s -1 ~| 

— 47r53(r,,) + — (5^^-3^^) {a.xp.yj (44) 



We have encountered here for the first time singular electron-electron operators. One can 
make them meaningful by appropriate regularization of the photon propagator, or by di- 
mensional regularization. In general it is a difficult problem and, as we have written in the 
introduction, the explicit solution was demonstrated for the positronium and helium atom 
only. 

6H3 is an correction that comes from 5*"^ term in Eq. (j22p . 



(a ■ pa ■ E + a ■ E a ■ pj . (45) 



16 

To calculate it, we have to return to the original expression for one-photon exchange. We 
assume that particle a interacts by this term, while particle b by nonrelativistic coupling 
e and obtain 

d^k 1 1 



m = EE 



a^b b 



(2 7r)4i p 16 m 



3 



Eq- Ho-k^ + ie 

+ (0|e^'^-^^^ — — (aa -Pa^a-k e'^-'^'^ + d^-""^ aa -Pa^a- k) 10) (46) 

Eq - Hq - k^' + I e ^ ' ) 
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We replace k —>■ —k in the second term, then perform the integral, and obtain 



' 2 i (2 7r)3^2l6m3 



■ Pa Oa " A^C^'^'^^" + C*^'^^'" ■ Pa ' k) {H, - E^) e-^^'^^l^) 

+ (0|e^^^-^^' (i/o - E,) (aa -Pa^a- ^e^'"^" + e*'"'" a, ■ Pa a, ■ Q 10)} (47) 

After commuting [Hq — Eq) with e^*'^'^'' one expresses this correction in terms of an effective 
operator 



^^3 = EE 



1 



aj^b b 



32 





■ 2 O! ' 




pI, 


Pa, 






I rabi 





(48) 



6H4 is the relativistic correction to transverse photon exchange. The first electron is 
coupled to A by the nonrelativistic term 



p.A-—a-B 

m 2m 



(49) 



and the second one by the relativistic correction, the 3'"'^ term in Eq. (|22j) 



(n^ — ea ■ B n'^ — n'^ ea ■ B) 



Q m-^ 

It is sufficient to calculate it in the nonretardation approximation 

a 



{p^2p - A + 2p - Ap^ + a ■ B a ■ B) (50) 



5^4 = E E 

ajtb b 

\pi , 1 



— + 7; — (o^b X Vfo 



'^PaP\+Pa (^a >< 
1 



2rafe 



r 



(51) 



-m 2 m 

It is convenient at this point to introduce a notation for the vector potential at the position 
of particle a which is produced by other particles 



« , ^ibrib\ Pb , a i^bX Tab) 



+ 



m ' 2m r^j, 



(52) 



then this correction can be written as 



5Hi = E UplPa ■ Aa + 2pa ■ AaPl + Pa " Va X A + " Va X A Pal (53) 



)m^ 



Let us notice that in the nonretardation approximation any correction can be simply ob- 
tained by replacing the magnetic field A by a static field Aa- We will use this fact in further 
calculations. 
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SH^ comes from the coupling 



8 m' 



a - {E X A-Ax E) 



(54) 



which is present in 4*^^ term in Eq. (j22I)- The resulting correction is obtained by replacing 
the fields E and A by the static fields produced by other electrons 

o2 



^ ^ 

SH.^ = V Ba- £a^ Aa- Aa^ £a 



5Hq comes from the coupling 



2 m 



■A' 



(55) 



(56) 



which is present in the second term of Eq. (j22I)- Again, in the nonretardation approximation 
the Aa field is being replaced by the static fields produced by other electrons 



^^6 = E^^1 



(57) 



5H-J is a retardation correction in the nonrelativistic single transverse photon exchange. 
To calculate it, we have to return to the general one-photon exchange expression, Eq. (|^. 
and take the transverse part of the photon propagator 



5E = -e 



1 



5'^ - 



fa{k) e^'^- ^ ^ ,n , , 3i{-k) e'^^^-'- 



(58) 



Eq- Ho - k^ + ie 

We assume that the product j^(fc) ji{—k) contains at most a single power of This allows 
one to perform the integration by encircling the only pole k^ = \k\ on '^{k^) > complex 
half plane and obtain 



5E 



d^k 



(271)^ 2k 



k' y 



1 



■3h{-k) e" 



-i k-ri, 



{a ^ b) . 



(59) 



Eo-Ho-k 

where k = \k\. By using the nonrelativistic form of / and taking the third term in the 
retardation expansion. 



1 1 Hn — En (Hq — EqY 



Eo-Ho-k 



k 



~y 



y 



(60) 



13 



where the first one contributes to the Breit-Pauli Hamiltonian, the second term to E^^\ and 
the third term gives 6Ei 



SEr 



EE 

a^b b 



1 



5'^ 



m 2m 



(61) 



,171 2m 

This is the most comphcated term in the evaluation, and we have to spht it into three parts 
with no spin, single spin and double spin terms 



5E7 = 5Ea + 5Eb + 5E, 



c 



(62) 



The part with double spin operators is 
5Ec 



EE-e^ 

a a^b 



(2 7r)3 2A;4 
One uses the commutation identity 



d^k (cxa X k) ■ (cxb X A;) / 



4m^ 



e*fe-. (i7g_^Q)2e-*fc-r, 0\ (g3) 



( k-ra 



2m^ 

to express this correction in terms of the effective operator 6Hc- 

a 



{Ho-EoY,e-^'^-''^ 



5H, 



c 



EE 



a>b 
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The part with no spin operator is 

^Ea = EE-e^ 



ay^b b 



dPk 
{271)^ 2k* 



2 . . 1 fri,Ti, 5^^ 



"irab 



k' y 



2 Tab V r^^ 



(64) 



(65) 



Pa(ik-r^ (ifo - Eo? e"*^-^" - (Ho - Eof} 
m ^ J m 

We subtracted here the term with k = 0. We ought to perform this in Eq. (jF)T|) . where lower 
order terms were subtracted, but for simplicity of writing we have not done it until now. 
We use another commutator identity 



N2l Pb 



(66) 



{Ho-EoYe-''-^^''-{Ho-Eo) 



{Ho - Eo) (e^^-'^- - 1) {Ho - Eo) + {Ho - Eo) 

n2 



Pb J h 



2 m 



■Tab 



- 1 



+ 



Pa 

2m. 



{Ho - Eo) + 



pI 

2m' 



pI 

2 m. 



(67) 



14 



and the integration formula 



to obtain the effective operator 6Ha 

a 



(68) 



SHa = EE 

a>b b 



b^,v-] """""" "' [VM] 



Pb '^ab'^ab - 3 



2 n 

ab 



2 m 



Tab 



'ribrib-^S'^r'i,, pi 



Tab 



2 m 



[V,pi] 



+P\ 



pI \r'abrib-^5'^ rl^, pi 



2 m 



Tab 



2m. 



Pb 



(69) 



The part with the single spin operator is 

[e'^^-'^ (Ho - Eof e-^^^-^^" aa x k ■ p, - p^ ■ ^b x ^ e^'"'" (i^o - i?o)' e"*^^'^^''} (70) 

With the help of the commutator in Eq. (jUTjl identity and the integral 

Airk 



d^k- 



A k-f 



k^ 



1 r 

2 r 



(71) 



one obtains 

5Hb 



EE 

a>b b 
-[Pa,V] 



a 



Aim? 



rab pI ' 



Tab 2 m 



[VJb] + 



pI 

2m' 



rab pI 



r-afe 2 m 



Pb 



2 '"ab 




[ pI 




X — 


-Pa 










.2m' 





rafe Pt 



Tab 2 m, 



Finally, the operator 5H-j is a sum of already derived parts 



(72) 



(73) 



is a retardation correction in a single transverse photon exchange, where one vertex 
is nonrelativistic, Eq. and the second comes from the 3rd term in Eq. 

e 



8 m' 



a ■ {E X p — px E) 



(74) 



With the help of Eq. one obtains the following expression for 5E^ 



SE. 



k' y 



k"^ J 16 m'^ 



I (e^ PaXBa+PaX Ba e''''''^ 



(75) 
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In the expansion of l/(-E'o "~ Hq — k) in Eq. ()60p the first term vanishes because h.c. and 
the second term is a correction of order ma^. After commuting {Hq — Eq) on the left one 
obtains the effective operator 5H^ 



ajtb b 



(2^ 



k'y\ 1 



k"^ J k"^ 16 



2m 



■ CTr, 



I e 
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•Aa ■ Pa ^ (^a + Pa ^ (^a ■ A 

a 1 P, 



6Hg is a one- and two-loop radiative correction 



- <?b X A;) + h.c. 



(76) 



(77) 



and its derivation requires a separate treatment. We base our treatment here on known 
results for helium, which in turn are based on hydrogen and positronium, and extend it to 
an arbitrary atom, as long as nonrelativistic expansion makes sense. 



where 



H 



Rl 



5:f!^ 9.61837 ^=(rj 



+ E E ^ (14.33134 - 3.42651 a„ • a^) S\rab) 



a>b b 
„2 



H 



R2 



0.17155 



EE — (-0.66526 + 0.08633 a, • a,) S^r^b) 



a>b 



m 



Z a 



Ara^ 



■ cr. 



+ EE 



a>b b 4 rl, 



<5^^'-3%^) [2a(2) + (a«)2] 

^ab ^ 



a 



4 r^f^ 



i^a + ^b) ■ Tab X {pa - Pb) 2 af^ 



(78) 



(79) 



« = 0.5, 

= -0.328 478 965... 



(80) 
(81) 
(82) 
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and K is the electron magnetic moment anomaly. 



IV. SUMMARY 



The obtained complete m contribution H'^^'> Eq. ()38|) is i n ag reement with the former 



derivation for the particular case of helium S 



m and P levels 



12| |. but is much more com- 



pact. Due to differing ways of representing various complicated operators, this comparison 
is rather nontrivial, and we had to refer in many case to momentum representation to find 
agreement. Since the present derivation differs from the former one, this agreement may be 
regarded as a justification of this, and as well as the former results. We have not derived here 
the term A in Eq. (jSj). It is obtained by matching the forward scattering amplitude in full 
QED with the one obtained from the effective Hamiltonian or NRQED, and it accounts for 
the contribution with high electron momentum. A depends however, on the regularization 
scheme, and once it is fixed it can be obtained in a similar way as in dimensional [2^ or 
photon propagator regularizations j2o| . 

can be used for high precision calculations of energy levels of few electron atoms, 
provided two difficulties are overcome. The first one is the algebraic elimination of electron- 
electron singularities. The elimination of electron-nucleus singularities was demonstrated 
on hydrogen and helium examples, and could easily be extended to an arbitrary atom. The 
elimination of electron-electron singularities using the dimensional regularization scheme 
was performed for the ground state of helium atom in 13], however this derivation was 
very complicated and so far this result has not been confirmed. The extension to more 
than two-electron atoms is even more complicated, therefore a new idea which will lead to 
elimination of electron-electron singularities is needed. The second difficulty is the lack of 
analytical values for integrals with basis sets, which fulfill the cusp conditions. For example, 



ipie 



for the Hylleraas basis set only three-electron integrals are known analytically [26|, |27 |. 
This cusp condition is necessary, because effective operators present in 5H^^^ contain many 
derivatives. For these reasons the calculation of ma^ contribution to atomic energy levels 
Kas been ^ccnpU^ed on,, fo. a few .ate. of he„u™ a... H Q H Q 



and hyperfine 



splitting 26 | , where a powerful random exponential basis set has been applied. For the 
three- 0, and four-electron systems the leading QED effects, namely the correction 
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of order ma^ have only recently been calculated. Due to developments in the Hylleraas 
basis sets we think the calculation of m contribution to lithium energy levels is 

now possible, particularly interesting is the Qia"^) correction to hyperfine splitting, which 
may be regarded as a benchmark calculations for MCDF or MBPT methods. Another 
interesting example is the fine structure of Pj levels, where electron-electron singularities 
are not present, such as for helium fine structure. 
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